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ABSTRACT 


The temperature stress day (TSD), defined as the difference in tem- 


perature of a well-watered and a water-stressed canopy, has been widely 


accepted. However, its dependence on vapor pressure (VPD) has never before 
been reported. The VPD dependence was proven both theoretically and 
empirically using the Penman-Monteith resistance equation and the non- 
water-stressed baseline equation relating canopy and air temperatures to 
VPD. The VPD effect on the canopy temperature difference was accounted for 
in a new plant water status indicator based on the TSD called the tem- 


perature stress day index (TSDI). The TSDI was defined as: 


where Т. is the measured canopy temperature, Тот, is the canopy temperature 
expected of a fully well-watered crop, and Tey is the canopy temperature 
expected of a fully water-stressed crop. Тот, and Tay were also defined 
theoretically and empirically and compared to measured Toy, and Tey values 
for alfalfa grown under both well-watered and water-stressed conditions. 
Calculated values closely approximated measured values indicating that the 
linear ratio of canopy temperatures used in the TSDI provides a good method 


for quantification of plant water stress. 








INTRODUCTION 

Detection of water ЕЕ by remote sensing techniques has been the 
subject of many studies beginning with Monteith and Szeicz (1962) and 
Tanner (1963). These scientists applied the then new science of infrared 
thermometry to the measurement of plant temperature which was reported to 
be an indicator of plant water status in the early 20th century (Miller and 
Saunders, 1923). Development of portable infrared thermometers since that 
time has led to the proliferation of their use. Subsequently, crop water 
stress indicies which utilize canopy temperature have been developed. Some 
of the indicies now in use include: the crop water stress index or CWSI 
developed empirically by Idso et al. (1981) and substantiated theoretically 
by Jackson et al. (1981), and the temperature difference between a well- . 
watered canopy and a water stressed canopy called the temperature stress 
day (TSD; Gardner et al, 1981). Operational use of the empirical CWSI 
requires concurrent measurements of canopy and air temperature and of the 
water vapor pressure deficit of the air. Use of the TSD requires con- 
current measurements of the canopy temperature of a well-watered field and 
of the field of interest. The TSD was said to be advantageous over CWSI in 
stress detection since TSD was thought to be independent of the environmen- 
tal factors upon which CWSI is based (Jackson, 1982). Still, the CWSI has 
become more widely accepted since well-watered fields are not always 
available as required by the TSD. Furthermore, there was no specified 
upper limit to which a TSD value might rise. Hence, no absolute level of 


stress could be calculated through the TSD. 
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This treatise will discuss the development of an upper limit for the 
TSD, and a crop water status index will be developed from those ideas. We 
will also show that the TSD method is, contrary to current opinion, depen- 
dent on the same environmental variables as the CWSI. The approach will be 
both empirical and theoretical. 
Theoretical Considerations 


The energy balance of a crop can be written 
R, = LE + H + G | 1 


where R, is net radiation, LE and H are the latent and sensible heat 
fluxes, respectively, and G is the soil heat flux (all terms expressed as 
W m2). Following the electrical resistance analog of Monteith (1973) and 


) modified by Verma and Barfield (1979), LE and H can be expressed as 


* 
pCp (ec - ea) 


LE = 一 -一 一 2 
Y (rc + Tay) 
(То-Та) 
0р as 7 


where p is the density of air (kg m3), Cp is the specific heat of air at 





constant pressure (J кеті С), y is the psychrometric constant (kPa ολ, 
ec is the seer ated vapor pressure (kPa) evaluated at canopy temperature 
(Te, units of C), eg is the actual vapor pressure of the air (kPa) eval- 
uated at the air temperature (Tg, units of C), гс is the canopy resistance 


(s m7!) and ray and ras are the aerodynamic resistances to the fluxes of 


water vapor and heat (s m 1), respectively. 














a4 
By substituting the identities of H and LE into Eqn. 1, setting G-0, 


and assuming equality of гау and ras, Tc can be determined by 


* | 
Rn Ya У(1+ге/та) (еа-еа) 


Tc = Та + p Cp "UNMNORr'rQG) — MN rU ra) 


It is obvious from Eqn. 4 that Т. depends on air temperature, net radi- 
ation, water vapor pressure deficit (VPD) of the air, and the aerodynamic 
and canopy resistances. 

From Eqn. 4, we can define an upper and a lower limit of Tc for 5 
given Та, Еп» ea, Ta and гс. For the maximum upper limit of Т. (Tey) at 
which no transpiration occurs, гс approaches infinity. Eqn. 4 then 


simplifies to 


Rn Га 
pCp 


Tcu = Ta + 
indicating energy exchanges of daly radiation and sensible heat. The 
result (Fig. 1) is a straight line relationship between Toy and Ta at а 
given Еп and га. For a crop well supplied with water, the potential lower 
limit of Тс (TeL) can be determined in а similar manner. This is done by 
setting τς to zero; the case of the plant acting as a free water surface. 


Eqn. 4 can then be expressed as 


| * 
En Ta Y (ea-ea) : 
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A family of curves can be generated from Eqn. 6 as a function of Tg and 
relative humidity (ед/е%). These curvilinear relationships are shown in 
Fig. 1 with relative humidity ranging from 0 to 100%. The most accurate 
solution to Eqn. 6 requires iteration by evaluating A at (TeL + Ta)/2. 
At 100% relative humidity, Τετ, approaches Τρι but remains 2 to 4 C lower 
than Toy. This occurs because canopy temperature becomes warmer than air 
temperature. Since the canopy is warmer than air temperature, the 
Saturated vapor pressure at the crop surface is larger than the saturated 
vapor pressure evaluated at То. A vapor pressure gradient still remains 
between the crop and the air. As a result, transpiration still continues 
in saturated air to cool the crop below Tey but not below Ta. 

The TSD can now be defined in terms of the controlling parameters 


which specify Т. and TeL, namely by combining Eqns. 4 and 6 such that 


3 (1+гс/та) 
TSD = T - Top = BIB М а "a 


* 1 1 
+ ες ps s 
(eaea) [с ЕЙТЕ | d 
When Т. is measured, Eqn. 7 сап be simplified to 


* 
TSD = T. - Tey = Т.------%--------- | 8 
c€ eh рб ` Aty AY | 
Since both R, and VPD vary diurnally and day-to-day, an identical value of 


TSD obtained from different days or even different hours may not 


necessarily indicate the same level of stress, For example, by altering 








relative humidity by only 25% (from 25 to 50%) at Ta=30 C, the lower 
canopy temperature limit rises from 18.5 C to 23.4 C--a difference of 
4.9 C. As a result, TSD changes by 4.9 C. This illustrates that TSD 
cannot be considered to be independent of environmental conditions. In 
order to be useful, TSD must be normalized in such a manner that these 
conditions can be accounted for. This process is described in the 
succeeding paragraphs. 

A maximum canopy temperature difference expected of a water stressed 
crop above that of a well-watered crop is determined from the appropriate 


combination of Eqns. 5 and 6 such that 





* 
r ‚ (еа-еа) 
cu cL oC AY AY 


Air temperature no longer appears explicitly in the formula, but appears 
implicitly in A e; and ед. The resulting solution, plotted now as a 
function of relative humidity, is shown in Fig. 2. The area under each 
curve is the range where measured canopy temperature differences are 
expected to occur at a specific relative humidity or VPD. 

The values obtained from Eqn. 9 represent the maximum theoretical 
canopy temperature difference which would be expected when a fully well- 
watered canopy is compared to a fully water-stressed canopy. However, the 
assumption that canopy resistance is negligible under well-watered con- 
ditions is rarely accurate (van Bavel and Ehrler, 1968). A more realistic 


representation is found by 








Several variables are required as inputs to Eqn. 11, some of which are 
difficult to obtain. However TcL is, in practice, concurrently measured 


with To. When this is done, Eqn. 9 simplifies to 


Te = TeL 
Кг 

п Га 
Та + T - TeL 


TSDI = 12 








In this simplified form, the TSDI requires only net radiation, air tem- 
perature and aerodynamic resistance values as well as the two canopy tem- 
peratures Т. and Το]. Тһе need for measurements of relative humidity or 
VPD and τρ], are obviated. Measurement of Ry, Та, Тс and Тот, are straight- 
forward while ra can be estimated by several methods such as those 
described in Jackson et al (1981) and Heilman and Kanemasu (1976) among 
others. 

The operational use of the theoretically derived TSDI can be 
illustrated using both the theoretical upper and lower canopy Papapa p | 
limits. By specifying Та, Rn; ra, and гс, a family of TSDI curves can be 
generated as shown in Fig. 3. For example, if Т.-27.0 C and VPD-3 kPa, 
thén' the TSDI would be 0.48. 

exque ч Empirical ΠΝ 

The empirical development of the TSDI follows closely the empirical 
development of the CWSI. The basis of the CWSI lies in the concept of a 
Е relationship of the canopy-air temperature difference (Т.-Та) as а 
function of vapor pressure deficit when the crop has ample water (Ehrler, 


1973). The basic form of the relationship is 





TeL - Ta = b+ ш(ез-е4) | 13 


where b апа п аге coefficients for the intercept and slope, respectively. 
All well-watered plant species studied to date have been — MM to ЕЕЕ 
bit a unique TcL-Ta vs. VPD relationship (Idso, 1982). The empirical CWSI 
was born from a ТЕ Е of these so-called "non-water-stressed baseline" 
relationships and the limit which Те-Та approaches in the complete absence 
of transpirational cooling (Idso et al, 1981). The upper limit, according 
to Pinter (1983) is a function of the well-watered baseline and can be 


expressed by 


Тол - Ta = b + ш(еа-еа+ь) 14 


where es rb is the value of the saturation vapor pressure at Ta + b. The 
lower and upper limit equations can be easily rearranged making direct 

calculation of Теа and TeL possible. When Eqns. 13 and 14 are combined, 
the maximum temperature difference between a fully water stressed Кен 


and a fully well-watered canopy can be determined such that 
Tou - τει, = m [(ел-еа+ъ) ux (es-ea)] 15 


This TSD, calculated empirically, also emphasizes the dependence of the 
‘canopy temperature difference on VPD. 
From these empirical relationships, the TSDI can be readily 
constructed in the following manner: 
* 
Te = Του les Ta =D + m(es-ea) 


TSDI = -人 _ = 一 ——a 16 
Teu - TeL m [(ea-ea+b) - (ea-ea)] 


) 
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An a priori knowledge of the non-water-stressed baseline is required before 
the TSDI can be empirically ай, In this form, the empirically 
derived TSDI and the empirically derived CWSI are identical. When Tc and 
Тор are Е only knowledge of the upper limit is required. 


Eqn. 16 then simplifies to 


T. - T 
TSDI - тысы > AR 17 
Ta 一 TeL + b + m(ea~eatb) | 


An example of the operational use of the empirical TSDI is illustrated in 
Fig. 4 for alfalfa (Meticago sativa L.) at Та-30 C. At а VPD of 3 ЕРА, а 
measured canopy temperature of 27 C would equate to a TSDI value of 0.37. 
MATERIALS AND METHODS 

Experimental verification of the TSDI was conducted from 15 Nov 1984 
through 09 Jul 1985 on a one-year-old stand of alfalfa (Meticago sativa L., 
су. Lew) grown in Phoenix, Arizona. The soil is a fine, loamy, mixed, 
calcareous, hyperthermic Anthropic Torrifluvent. Тһе 0.5 ha field was 
divided into 18 plots ranging in size from 160 to 325 m2. The plots were 
arranged in 6 blocks of 3 plots each and 3 irrigation treatments were ran- 
domly imposed on two of the blocks during each growth cycle. The 5 growth 
cycles were determined by the cutting dates, which were 15 Nov 1984 and 12 ' 
Feb, 28 Mar, 07 May, 06 Jun, and 09 Jul 1985. The irrigation treatments 
were designated Early, Late and Dry and were designed to distribute the 
imposed water Seress СТЕ: different periods of a growth cycle. The 
Early irrigation treatment was irrigated within 3 days after harvest to 


field capacity at a depth of 0.9 m, but received no further supplemental 
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irrigations. A shallow irrigation was an assurance of stress development 
later in the growing cycle. The Late treatment was irrigated to field 
capacity at a depth of 3 m halfway through the growth cycle but received no 
initial irrigation following harvest. The Dry treatment received no irri- 
gation water from harvest to harvest. Following harvest, the irrigation 
treatments were rotated to two new blocks to avoid any cumulative stress 
affects on stand density. In addition to the three deficit irrigation 
treatments, two Wet or fully irigated treatments were maintained in two of 
the remaining non-treatment blocks. These plots were irrigated immediately 
after harvest and whenever the soil water content, as determined by thrice- 
weekly neutron probe measurements, declined to 507 of field capacity in the 
top 3 m of soil. 

Canopy temperature (Everest intérsctenes i. model 110, 4? FOV) and 
concurrent air temperature, vapor pressure deficit (Bendix ο στ, 
model 566) and net radiation (Fritschen; 1963, 1965) measurements were 
made daily at 1330 hrs MST under full sun and ground cover conditions. 
Canopy temperature of each plot was measured from the four cardinal 
directions was determined at six different sites per view direction. The 
infrared thermometer was held at a declination angle of 30° from the hori- 
zontal resulting in a spot size of about 0.25 m2. Psychrometric obser- 
vations collected both immediately preceding and following the canopy 
temperature — were averaged to give mean air temperature and 


vapor pressure deficit. Net radiometers placed in 6 plots (both treatment 
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and nontreatment plots) at a height of 1.5 m were scanned automatically 
every 1.5 min during canopy Еее measurements. These data were time 
averaged over a 30 minute period to give a mean net radiation value. 

Study-long empirical and theoretical values of Тот, and Teu were calcu- 
lated and compared to measured values. Non-water-stressed baselines for 
each growing cycle and for the entire study were computed according to 
Eqn. 13 from the daily midday observations of canopy and air temperature 
and vapor pressure in the Wet plots. Im addition, a single diurnal obser- 
vation of these variables, recorded every 30 min from 1000 through 1600 
hours MST, was also included in the computation of the non-water-stressed 
baseline for each growth cycle. These baselines were compared to those 
reported by Idso et al (1981) for alfalfa. 

RESULTS AND DISCUSSION 

The six month period of measurement provided an opportunity for large | 
ranges of air and canopy temperature, net radiation and vapor pressure 
deficit to be observed through five growth cycles. Air temperature at the 
time of canopy temperature measurements ranged from a minimum of 14.7 to a 
maximum of 40.3 C during the course of the study. Concurrent values of net 
radiation ranged from 310 to 720 W ш-2 and VPD ranged from 0.89 to 6.59 
kPa. The response of measured canopy temperature to these three environ- 
mental parameters and to the imposed water stress is evident in Fig. 5. 
These data are plotted for only clear sky conditions and when canopy cover 
was 100% developed. Maximum values of Tc at a given Ta represent Toy while 
minimum values of То represent Тор. Little spread was observed in canopy 


temperature at air temperatures less than 25 C. These air temperatures 
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represent the winter months when precipitation events occurred frequently 
and no water stress developed. Hence, no Тс. values were recorded during 
this time period. The wet climate and cool air temperatures were also con- 
ducive to the growth of downy mildew (Peronospora trifoliorum d By.) and 
blue alfalfa aphids (Acyrthosiphon kondoi Shinji) and these two pests were 
prevalent during the Jan-Feb growing period. The probable result was a 
warmer Wet canopy temperature than would otherwise have been expected. 
Therefore, no values of TeL or Teu were recorded during the winter months. 
No infestation of either mildew or aphids was observed during the warmer 
and drier period of Mar-July.' Rainfall during the latter period (T. > 
25 C) was also insufficient for optimum growth and severe water stress was 
visually inferred from plant wilting and leaf abscission. The large spread 
in canopy temperature was recorded at air temperatures >25 C and differen- 
ces (Dry-Wet) as large as 10 C were measured. Maximum and minimum T, 
during this time period were considered to be representative of Tcu and 
TeL, respectively. 

Empirically derived TcL and Tcu are shown for the entire study in 

Fig. 5. In order to plot Тот, and Toy ав a function of Ta, VPD was first 
ee as a function of Ta as shown in Fig 6. A simple second order 
polynomial equation with an R2 of 0.980 adequately specified the T,:VPD 
relationship. With the T4:VPD prediction equation, Tey, was computed using 
Eqn. 13 with the slope and intercept values for alfalfa reported by Idso et 
al (1981). Empirical Teu was obtained in a similar manner using Eqn. 14. 
The calculated values of Тот, and Toy fit closely the measured minimum and 


maximum canopy temperatures. However, empirical estimates of Тот, and Tey 
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were in closer agreement to measured values when slope and intercept 
coefficients, calculated from Eqn. 13 with Wet data, were used. These 
coefficients are summarized in Table 1 for each growth cycle and for the 
entire study. Тһе study-long ЕВЕ verë also used to calculate 
empirical Тот, and Toy, which are plotted in Fig. 5 along with the TeL and | 
Toy derived from Idso et al (1981). | 

The general seasonal trend ёте in the coefficients was for the 
slope to become decreasingly negative and for the intercept to decline, 
even becoming negative as the season progressed. This finding is in 
opposition to Idso (1982) who stated that no negative intercepts had been 
recorded for 26 different plant species, among them alfalfa. The slopes 
and intercepts obtained during the two growth cycles between 28 Mar and 06 
Jun most closely approximated the Idso et al (1981) baseline. It was 
during this same time of year that the bata for their baseline, grown at 
the same location, was collected. | 

The time-of-year dependent o shown for this study were well 
correlated to VPD (see Table 1). The changing slopes and intercepts 
suggest а curvilinear relationship of T¿-T¿ with VPD. This is precisely 
the finding of Jackson et al (1981) in their, theoretical approach to the 
calculation of Т.-Та. Kirkham et al (1983) also indicated that the range 
of VPD greatly affected the slope and intercept of the non-water-stressed 
baseline. Indeed, regression analysis of a second order polynomial fit to 
our data showed the higher order equation to be significant (P<0.05) with 
the equation being Т.-Та = 4.04 - 3.17(VPD) + 0.142(VPD)2. Coefficient of 


determination for the fit was 0.937. 
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We conclude, therefore, that different coefficients for the linear 
relationship of Т.-Та vs. VPD сап be obtained in different seasons or in 
different climates, depending on the range of VPD over which the relation- 
ship is derived. Ме conclude further, that net radiation plays only a 
minor role in the Т.-Та vs. VPD relationship as demonstrated by the good- 
ness of fit of the calculated empirical Ter. shown in Fig. 5. m finally, 
the empirical approach to the calculation of Το], and Ίσα is sufficiently 
unease to justify its use. | 

Once TeL and Teu are determined empirically, they can be used in the 
TSDI. Any measured Т. will fall between hese two temperature ng 
The TSDI, as defined by Eqn. 16, assigns a relative stress level according 
to the distance Т. is removed from the Тот, baseline. 

Theoretical Тот, and Toy, plotted together with даба e T the study, 
are shown in Fig. 7. Toy was obtained from Eqn. 5 by first setting а 
nominal value for Rn at 600 W m~2 and then iterating to find ra such that 
_the line closely approximated the maximum observed value of Те at a speci- 
fic Та. A similar approach was described by Jackson et al (1981) who 
applied the technique to determine rg in wheat (Triticum durum Desf.). The 
ra value obtained for alfalfa using this approach was 6s ml, While this 
value appears to be small, it is not poussent with ra values for other 
crops which were obtained in a similar manner. Jackson et al (1981) 
reported a value of 10 s ml for wheat and Idso et al (1986) reported а 
value of 24 s ml for water hyacinth (Eichhornia crassipes (Mart). Solms). 


Since ra is largely a function of canopy characteristics, and with leaf 


size being a primary factor of any canopy characteristic, it follows that 
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the lowest value of ra would be reported for alfalfa and that water 
hyacinth would have the largest га. Wheat, with an intermediate leaf 
size, would fall somewhere between the two extremes. 

After computation of τη, гс was determined іп a similar manner to 
that described for rg using Eqn. 4, The Ta:VPD do previously 
described was used to obtain values of VPD at a given T, and a уййне of 
600 W m"2 was used for net radiation. Canopy resistance was determined 
iteratively by substituting values of r< into the equation to make the line | 
fall through the minimum То values, A value of 15 s m | for г. provided 
the best fit. The iterative calculation of rg and r< provided a good fit 
of TeL to the data for the entire range of Tg. Тһе same was true for 
Toy during the time period when stress was evident (T4525 C). 

Net radiation enters into each equation for TeL and Toye Since Rn 
was not constant throughout the study, the upper and lower values of Rn 
(310 and 720 W m2, respectively) were used to indicate the extremes of 
Toy and To]. Changing net radiation by approximately 400 W m2 altered 
Тот, by only 1.0 C and Toy was altered by only 2.0 C at a Ta of 40 C. 

This indicates that net radiation served only a trivial role in the 
determination of Teg; and Teu» 

Once appropriate values of rg and r. к been determined, and the 
theoretical values of TeL and Toy have been calculated from measured 
values of Та, VPD, and Ry, a relative жене value сап be assigned to 


any measured Тс via the theoretical TSDI. 











17 
CONCLUS IONS 

The temperature stress day was proven, both empirically and theoreti- 
cally, to be dependent on VPD. As such the TSD is unacceptable for use in 
stress detection when a. constant VPD cannot be maintained. Instead, the 
TSDI was proposed as a means to account for this dependence. In order to 
calculate a TSDI, some knowledge of the expected minimum and maximum canopy 
temperature extremes at a given air temperature and vapor pressure deficit 
is required. Both theoretical and empirical methods for the derivation of 
these two extremes were presented. The calculated extremes were compared 
to measured extremes and both the theoretical and empirical values were in 
good agreement. As such, the TSDI appears to be a viable approach to 
remote monitoring of water stress. 

Some limitations to its use are evident for both the theoretical and ` 
empirical approaches. In the case of the theoretical approach, values of 
га and гс аге, in practice, difficult to obtain. Inaccurate estimates of 
either or both parameters alter the predicted TeL and Teu, thus affecting 
the TSDI. In the case of the empirical approach, an a priori knowledge of 
the baseline coefficients is needed before the TSDI can be used. However, 
many of these limitations can be avoided if the TSDI is applied to thermal 
data obtained either from either aircraft sti: Minimum canopy 
temperatures for the TSDI would not uud to be calculated in this case, but 
could be determined from measured Тс values. The coolest canopy tempera- 
tures measured from an airborne platform could be assumed to — n 
Тег» particularly in an irrigated area. This assumption would obviate 


the need for r< and VPD, the former being particularly troublesome to 
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determine. Only maximum canopy temperatures would need to be estimated 
according to Eqn. 5, and R, can be determined from the same platform 
(Jackson et al, 1985) as ie temperature data. Ground measurements of Tg 
could be obtained from automated weather data collection networks similar 
to that in operation in Nebraska (Hubbard at al, 1983) and estimates of ra 
ЕТ be obtained in a similar manner reported here. То, could then be 
rapidly calculated. Only this simplified data set would be needed to 
calculate TSDI, as specified in Eqn. 12, which could then be used to assess 
the relative water stress for entire fields or farms. Irrigation sche- 

_ duling on a grand scale based on canopy temperatures could be adopted. 
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Table 1. Results of linear regression analysis for T,j-T4 vs. VPD. Temperature is in units of Celcius and. VPD 


is in units of kPa. 
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Time period slope intercept R2* n Syex 
15 Noy I2 ве diseased canopy 
12 Feb - 28 Mar -2.64 4.06 0.763 47 0.557 
28 Mar - 07 May -1.87 1.43 0.888 102 0.606 
07 May - 06 Jun -1.99 1.01 0.883 42 0.436 
06 Jun - 09 Jul -1.60 -0.33 0.730 81 0.805 

All data -2.02 2.04 0.926 272 0.851 


0 


х R2= 


n 


Syex = 


51 


coefficient of determination 

number of observations 

standard error of the estimate of y on x 
standard error of the slope estimate 


standard error of the intercept estimate. 


0.220 
0.066 
0.114 
0.109 


0.035 


VPD range 

1.26 2.47 
1.66 - 4.84 
2.27 4.65 
3.69 6.59 
1.26 6.59 














Fig. l. 


Fig. 2. 


Fig. 3. 


Fig. 4. 


Fig. 5. 


Fig. 6. 


Fig. 7. 


LIST OF FIGURES 
Theoretically derived upper and lower canopy temperatures 
specified by Eqns. 5 and 6 where R, = 600 W m^2 and Ta = 
10 s wl, 
Theoretically derived canopy temperature difference (Те.-Тст,) 
specified by Eqn. 9 where R, = 600 W m~2 and ra = 10 8 ml. 
Theoretically derived TSDI specified by Eqn. 11 where R, = 600 
Wm 2, ra 10 s ml, το = 0 8 ml and T4 = 30 С. 
Empirically derived TSDI specified by Eqn. 16 using coefficients 
from the non-water-stressed baseline reported by Idso et al (1981) 
at T4 = 30 C. 
Alfalfa canopy temperatures obtained under clear skies at 1007 
ground cover from all irrigation treatments for the period 15 Nov 
1984 through 09 July 1985. Upper and lower limit lines denote 
empirically derived Teu and TeL from the non-water-stressed 
baseline reported by Idso et al (1981) and from the non-water- 
stressed baseline determined for these data. 
Vapor pressure deficit dependence on air temperature as specified 
by a second order polynomial equation. . | | 
As in Fig. 5 except upper and lower limit lines indicate 
theoretical Teu and Тот, determined from iteration of Eqns. 4 and 
5 by setting Rn = 600 W m”2 апа usine the T4:VPD dependence as 


shown in Fig. 6. 
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COMPARATIVE EVALUATION OF 
THE IDSO, JACKSON AND STANDARD OIL 
WATER STRESS INDICES 


By 
Bronson Gardner 


The purpose of this report is to summarize the theoretical behavior of 
three methods for calculating crop water stress index (CWSI) values. The 
basis for this analysis stems from field observations made in 1986 field 
trials which indicate that under some circumstances the current method ` 
of calculating CWSI (the Idso method) yields results which are difficult to 
interpret. Consequently, an intensive theoretical analysis of three 
methods of computing CWSI was developed. Two of the methods (one by 
Idso and one by Jackson) have been published in the literature. The third 
technique is a proprietary method developed at Standard Oil. 


Idso Index 


| .|dso et al. (1981) proposed a technique for evaluating the severity of 
water stress exhibited by a crop. The technique was based on the 

_ assumption of a linear relationship between the minimum сгор-аіг 

temperature difference (Delta-T) and the vapor pressure deficit of the air 


(VPD): 


| Min = a + D'VPD | (1) 
A crop water stress index (CWSI) was defined as 


CWSI = (Tc-Ta)-Min/(Max-Min) (2) 


where Tc = crop temperature 
Ta= air temperature 
max= theoretical maximum crop-air temperature difference 
Max was calculated by a linear extrapolation based of the Min line into the 
negative VPD region. 

















. This approach assumes that the coefficients a and b are independent | 


of net radiation (Rn) and air temperature effects. Evaluating the  Ὃ 5. 


coefficients for a given crop can be a difficult task , often requiring — ( 
several days worth of data, although there are times when meteorological 
conditions provide a sufficient range in VPD to determine a and b from 
hourly data collected on a single day. Іп spite of these difficulties, the 
Idso approach is appealing because of it's computational simplicity and 
because no windspeed or radiation terms are included in the index. 


Jackson Index 


By contrast, Jackson et al. (1981) proposed a CWSI based on the 
energy balance equation. Jackson et al. (1981) proposed that Min be 
calculated аз: | 


Min = га*Вп *а(1+гс/га) - УРО 
Rho*cp *d+g(1+rc/ra) d + g(1+rc/ra) 


where ra is the aerodynamic resistance of the crop , rc is the canopy 
resistance to vapor transport, g is the psychrometric constant, d is the 
slope of the saturated vapor pressure curve, VPD is the vapor pressure of 
the air, Rn is net radiation , Rho is the density of air, and cp is the heat 
capacity of air. Max was calculated as : | 


Max = ra*RN/(Rho*cp) (3) 


Jackson et al. (1981) recommended that the ratio rc/ra be evaluated with 
the following equation : 


rc/ra = [а*га*п/(Вһо*ср)] - (Tc-Ta)*(d+g)-VPD (4) 
ο" [(Тс-Та) - ra*Rn/(Rho*cp)] 


After rc/ra was calculated using this procedure, CWSI was calculated : 


CWSI = g * ( 1«rc/ra)- q*(1+rcp/r (5) 
d + g*(1 + rc/ra) 


where rcp is the crop resistance at potential transpiration rates for a 

given crop. lt can be seen that algorithms are also needed for estimating 
aerodynamic resistances and rcp values, in addition to the necessity of 
measuring net radiation and windspeeds. For this analysis, rc/ra was 
assumed to be 0.5 and rcp was set equal to O. 














- Standard ОЙ has developed a proprietary index which serves as an 
alternative to the Idso and Jackson equations. Evaluation of this index 
requires inputs of Tc, Ta, relative humidity and net radiation as a percent 
of the maximum expected net radiation. Evaluation of the coefficients in 
this index also requires a knowledge of the maximum crop-air temperature ` 
difference for a given crop at high air temperatures. 


Comparison Between the Indices 


In order to evaluate the overall performance of the Idso, Jackson and 
Standard Oil indices, it was necessary to generate a standard data set. 
In this standard data set, Rn was varied between 50 and 750 wim, in 
units of 50; Ta was varied between 15.6 C (60 F) and 43.3 C (110 F) in 
steps of 5.6 C (10 F); Tc-Ta ranged between -7.5 C to 22.5 C in steps of 
2.5 C ; VPD ranged between 0 and 5 kPa in steps of 0.2, providing the VPD 
value at a given Ta did not exceed the saturated vapor pressure, errors in 
measuring Tc (Ec) or decreases in Tc due to changes in Rn ranged between 
1 and 3 C. This data set contained 15210 combinations of Rn, Ta, (Tc-Ta), 
VPD and Ec. It should be remembered that the theoretical ranges in all 
three indices is between 0 and 1. 

The Idso index was evaluated for corn by using coefficients given in 
Idso (1982), with Tc-Ta limited to a maximum of 12 c. This standard data 
_ get showed that the Idso index is well behaved above VPD values of about 
2.0. Below VPD values of 2.0, the Idso scale rapdily expands to -12 to 
+12. The consistency of the scale above VPD's of 2.0 makes the 1050 
scale useable within these limitations. In practice, this would mean that 
the meteorological conditions under which the Idso index could be used 
accurately are well defined (Table 1). | 

The Idso index was also extremely sensitive to errors in measuring 
Tc or to changes in Tc as a result of fluctutions in Rn (figure 2). Тһе 
scale on this index also is a function of the magnitude of the value of 
Tc-Ta (figure 3). These results show that unless the Idso index is used in 
relatively hot dry climates (high VPD values) and unless Tc , Ta, and VPD 
are measured accurately under sunny conditions that significant difficulty 
in interpreting Idso's CWSI can occur. | 











Table 1 . Approximate regions where the Idso index provides consistent 


results. 


Air temperatures and 
must be greater than 


68.0 F 
75.2 F 
82.4 F 
89.6 F 
96.8 F 
104.0 F 
111.2 F 


Relative humidity ` 
must be less than 


14% 
33 96 
47 % 
58 % 
66 % 
73% 
78 % 

















The behavior of the Jackson index is similar to that of Idso (Fig. 4 -7). 
Jackson's index is a strong function of net radiation and vapor pressure | 
deficit. The scale of the index ranges between -8 and +14 at low VPD 
values. The index is an improvement over the Idso index, however. Note 
that the errors caused by decreases in canopy temperature (or errors in 


Tc) are more systematic than those of in the 1450 equation. 


The foregoing analysis is the basis for the new CWSI index . This new ` 
index has a nearly flat scale between over the entire net radiation and 
vapor pressure deficit ranges on which it was tested (Fig 8-9). As Tc 
changes due to errors or reductions in Rn, index values decrease bya 
maximum of only 0.4 ( as compared to several times 0.4 for the Idso and 
Jackson indices) (Fig 10). The relationship between the new CWSI and 
Tc-Ta is linear (Fig 11). At any given value of Tc-Ta, there is a clearly 
defined minimum and maximum value of CWSI. The exact value is then 
dependent on net radiation and VPD. Research using this new index has 
begun. 
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RESPONSE OF IDSO INDEX TO ERRORS/CHANGES IN CROP TEMPERATURE 2 қ 
СКОР ТЕМР ERRORS CORN OF 1 TO 3 C 
CROP = CORN 


SYMBOL IS VALUE OF ERROR/CHANGE FIGURE 2 
SYMBOL IS VALUE OF ERROR/CHANGE 
SYMBOL IS VALUE OF ERROR/CHANGE 

1:1 LINE -- SYMBOL USED IS » 





##111111111232 3 3 





+ 
| *###111111112221233222 
##111111121211323223332 
1.2 + *#111111111213222323 23 32 
| ###11111111222222232232 3 2 
»«11111111112122 232 3 2 3 2 
1.9 + ###1111112222222232 2323 3 3 ` 
| wee 2221211122232232 2 3 2 332 3 
«»1111111111113121232323 3 3 3 3 
0.8 + +#%11111111111332123331 3 2 3 3 
I | жєж 111111213112131 312 331 3 3 
D *+*111111112122132132 33 313 3 3 
S 9.6 + ** 111111111212 122 232 13 3 3 
0 | **##1111112112221322 22 1 3 3 31 
*«11111111211221213231232 23 1 3 
I 8.4 + **#*11111111212321232323232 2 13 3 
N *#*#1111111212212123223232331 3 23 
D %4111111112121 22 221 33132 232 32 
E 8.2 + ##%#11111111121212312 23232 33 3 
X »»111111111121223232 3 33333332 2 2 
| %%%11111111122222323232 2 322 2 2 
М 6.0 + *%+11111222222222723 3 3 3 33 3 3 3 3 3 3 3 
I *##%111112111112222 232 2 2 2 
T з%%11111112222333 322333 3 3 3 323 2 2 2 
H -8.2 + ###111211221113 222232323 з 333 32 2 2 
wee 11112112 2313 3332 33 2323 3 3 
E | | 中 中 中 11111222122 322 333 3 2 3 
R -0.4 + әже 111111123213 23 32 33 3 2 3 3 
R жеж 11111122133213 2 3 23 . 3 3 2 з 3 
0 | ++» 111111112123 3133 23 33 3 3 2 3 
R -6.6 + 中 中 中 11111121222123 32 3 223 3 3 3 2 3 
s www 11111122222132 32 3 2 3 3 | 3 2 
| wee 11121111212221 2 2 3 3 
-0.8 + eee 111111212212 213 2 2 3 33 3 3 2 
әзі 11111112223123221 23 2 3 3 ‚ 3 
| i a he ale 11111212312223 313 3 2 3 3 
-1.6 + eee 1111112112132313 3233 2 2 33 3 
Ii: 111111222323213223 33 23 2 3 3 
| жж 1111122222 22 3 323 3 3 3 
-1.2 + Ii. 111122212231232132 3 323 3 
жеж 11111211222232 23 32 23 3 
| II. 1111 11211233313 21 2 3 3 3 3 
+ 中 中 11 1 22 233 23232 3 3 2 2 
νον αι ως ------------ 十 一 一 一 一 一 一 一 一 一 一 一 二 一 一 一 一 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 ~ 中 一 一 一 一 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 一 二 一 一 一 一 一 一 一 一 一 一 一 +- 
-2.0 -1.8 -1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 1.6 2.0 


IDSO INDEX WITHOUT ERRORS 
2607 OBS HAD MISSING VALUES OR WERE OUT OF RANGE 3923 OBS HIDDEN 














RESPONSE OF IDSO INDEX TO TC-TA 


FIGURE 3 
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JACKSON INDEX VERSUS VPD 
CROP=CORN FIGURE 5 
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RESPONSE OF JACKSON INDEX TO A REDUCTION IN RN | 1 
RN DECREASES FROM 600 TO 200 W/M2, CROP TEMP DECREASES FROM 1 TO 3 C 
CROP=CORN 


SYMBOL IS VALUE OF ERROR/CHANGE IN TC FIGURE 6 
SYMBOL IS VALUE OF ERROR/CHANGE IN TC 
SYMBOL IS VALUE OF ERROR/CHANGE IN TC 

1:1 LINE-- SYMBOL USED IS s 


1.4 + 111 33 
| 111 33 
| 2222 111 
1.2 + 2222 
| 2222 
1.0 + 33333333 中 
1111111 
11111111 中 中 
J 0.8 + 2222 II: 
A 22222 3 жж 
С 2212 3333 * 
K 0.6 + 1111 3333 TI. 
S 1111 22333 中 中 中 由 
0. 1111122233«s»*« 
N 0.4 + 1 11122223308 
11 1111223335 
I 1111222333 
N 9.2 + 11111222333 
D x 11111222333 
E | 11111222333 
x 6.6 + 111112222333 
1111112222333 
w | *###111112223333 
I -6.2 + 中 中 中 111112223333 
T 中 中 中 中 1 1111122233333 
H | 中 中 中 由 11111112233333 
-0.4 + II. 11111111222333 
E III. 11111111222333 
R | ean 1 11111222233 
R -0.6 + II. 11111112232333 
0 жй 1 1111112222233 
R | 中 4 1 11111122222 3 
S - + * 1121111222 333 
1 11111112223 33 
| i 1 11111122333333 
-1.0 + 111111122233 33 
1 11111112233 3 
| 1 111211222333 
-1.2 * 1 11111322333 
21 111222233 
| | 211122223233 
-1.4 + 1 1112 3 23 
= =» «ке «е ою — A = === --+-- = === === ---+-------+------- - av «к «ә «м απο απο ———————— ofh o ——— ——Ó—Ó————— cu ds e l m al e 
-1.4 -1.2 -1.9 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


JACKSON INDEX WITHOUT ERRORS 
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RESPONSE OF NEW INDEX TO A REDUCTION IN RN 1 
RN DECREASES FROM 600 TO 200 W/M2, СВОР TEMP DECREASES FROM 1 ТО З С 
CROP=COR 


SYMBOL IS VALUE OF DECREASE FIGURE 10 
SYMBOL IS VALUE OF DECREASE 
SYMBOL IS VALUE OF DECREASE 
1:1 LINE -- SYMBOL USED IS » 








1.4 | 
1.2 + 
| 1 
1 
1.6 + 2 
| 3 
11111 3 
6.6 + 2222 
| 11 22233 
М 11112 333 
Е 6.6 + 11111222 3 
W | 111222233 
1111122233 
I 0.4 + 111112223333 
N | 111111222333 
D 11111222233 
E 9.2 + 111111222333 
x | 1111112223333 
111111122233 ` 
W 9.0 + 111111222233 
I | 1111111222333 
T 111111122223333 
H -ὐ.2 + 1111111122233 
| 11111111222333 
E 11111111223333 
R -6.4 + 1111111112233333 
R | 1##111111222333 
0 | 111 1112222333 
R -6.6 + 111111112223333 
s . | 1 111111112223333 
* 1 111112222333 
-0.8 + 11 2 2222222333 
| | 2 2 23333333 © 
i 2 33333 
-1.8 j 333 
-1.2 | 
-1.4 + 
一 一 一 一 一 一 +-+ 二 ~ 一 一 -4 一 一 一 一 = 一 -4 一 一 一 一 一 一 一 一 一 一 一 一 一 一 卡 一 一 一 一 一 一 一 直 一 一 一 一 一 一 一 二 一 一 一 一 一 一 一 和 一 一 一 一 一 一 一 个 一 一 一 一 一 一 一 个 一 一 一 一 一 一 一 上 一 一 
-1.4 -1.2 -1.9 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


NEW INDEX WITHOUT ERRORS 
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RESPONSE OF NEW INDEX TO TC-TA 


FIGURE 11 


CROP=CORN 
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12.5 


一 一 十 一 一 一 一 一 一 一 一 一 一 一 一 一 一 二 一 
10.0 


0.0 


-5 . 0 


| 
-7.5 


-2.5 


DIFF 
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^ fi B. R. Gardner Copy to: 


Improved Technique for Measuring Non-Stressed Baseline 
Values | 


I believe we should incorporate an alternate index into 
our options. There are users who complain about 
inconsistency in the current index. Inconsistencies are 
the result of errors in estimating the true minimum 
temperature of a crop for a given day. The factors which 
may cause these errors include: | 


Factor 


Sunlight variability 

Percent ground cover 

Windspeed variations 

CO5 variations | | 
(This factor, just recent discovered) CO» varies 
with soil moisture and the amount of decomposing 
organic matter. | 

Age of Crop | 

Sampling technique 

Time of day 

Genotype 

Surface wetness of the soil 

Instrument differences 

Air temperature fluctuations 

Baseline errors 


It is probable that no mathematical model will be 


developed in the near-term which can remove all these 


sources of variability. It is also clear that Scheduler 
sales/reputation hinge on the ability to absolutely 
minimize all sources of variability not associated with 
transpiration. In other words, the Scheduler index must 
be consistent despite all aforementioned potential sources 
of error. | | 


There is a way to accomplish this objective. I have 


tested this method on the cucumber data from Wisconsin and 


am convinced that this procedure should become a part of 
our options in selling/using the Scheduler. 
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The method requires a user to identify an area which has 
been recently irrigated (a piece of information which 
requires very little judgment). The user would simply 
walk into this area and take 30-40 samples. No index 
would be produced at this point. The user has just 
calibrated his baseline to adjust it to all the 
aforementioned sources of variability, including 
instrumentation. 


How is this done since the stress status of the field is 
not known (but by definition it should not be in a 
severely stressed status)? After the samples are taken, 
the CPU calculates the mean and standard deviation of the 
crop temperature and the mean sunlight level. The 
quantity [(Tc-Ta) - Std] for a reasonably non- stressed 
field is the minimum baseline. 


How would this change the way in which the current index 
is calculated? | 


The current index is: 
10 * (Diff - Min)/ (Max - Min) 
where Min = a + b * VPD 
The new index is: 
10 * (Diff - NMin) / (Max - Min) 
where NMin is the statistical minimum 
AS you can see, the denominators are identical, since the 
denominator represents the potential spread between 
Stressed and non-stressed plants. All that has been done 
is to substitute a statistically measured value for Min in 
pace of a computed value based on VPD. The suggested 


substitution can be expected to greatly improve the 
stability of the index. 
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